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Abstract 

There is considerable interest in determining the optimal polarizations that maxi- 
mize contrast between two scattering classes in polarimetric radar images. In this paper a 
systematic approach is presented for obtaining the optimal polarimetric matched filter, i.e.. 
that filter which produces maximum contrast between two scattering classes. The max- 
imization procedure involves solving an eigenvalue problem where the eigenvector corre- 
sponding to the maximum contrast ratio is optimal polarimetric matched filter. To exhibit 
the physical significance of this filter, it is transformed into its associated transmitting and 
receiving polarization states, written in terms of horizontal and vertical vector components. 
For the special case where the transmitting polarization is fixed, the receiving polariza- 


tion which maximizes the contrast 


ratio is also obtained. 


Polarimetric filtering is then 


applied to synthetic aperture radar images obtained from the Jet Propulsion Laboratory. 
It is shown, both numerically and through the use of radar imagery, that maximum image 


contrast can be realized when data is processed with the optimal polarimetric matched 


filter. 
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Contrast enhancement is a processing technique which modifies the input data struc- 
ture so that either the human observer, computer, or other hardware devices can extract 
certain information from the processed data more readily after t lie change Huany. 1969 . 
In this paper, the polarimetric properties of the radar return are utilized to enhance the 
contrast between two scattering classes. It is assumed that complete <i priori statisti- 
cal knowledge of the two scattering classes or types exists and the polarimetric signals 
backscattered from the two scattering classes are independent C'han, 1981 . The process- 
ing requirement is then to determine the optimal transmitting and receiving polarization 
state which will maximize the separation of the average power returns between the two 
classes. Applying such a technique to radar imagery will allow for better discrimination of 
the two classes. 

For two deterministic scatterers. completely characterized by 2x2 complex scatter- 
ing matrices. Kozlov 1979' introduced a method for computing the optimal polarization 
state which involves transformation of the scattering matrix of each of the two objects into 
a preferred polarization basis. Solutions involving scattering matrix co-pol and cross-pol 
(polarization) nulls have also been presented by h'cnnaugh 1949). Chan 1981’, Boerncr 
1982 . ct al.. Micros, 1983; Ncspor ct al., 1984; McCormick and Hendry, 1985]. More 
recently, Kostinski and Boerncr [1987] determined the transmitting and receiving polar- 
ization state which produced maximum contrast between two classes represented by their 
Graves power matrices. This technique involved maximizing the expected power return 
from one deterministic scattering class with respect to another. 

Although these are viable procedures when dealing with deterministic scattering 
classes, they can not be utilized in the case of statistically distributed scattering objects, 
e.g.. terrain clutter. A deterministic scatterer can be characterized bv a scattering matrix, 
whereas nondeterministic scatterers must be represented either by an average Mueller 
matrix or equivalently, by a polarimetric covariance matrix, which are the second order 
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statistics of the scattering matrix. To enhance the contrast between nondeterministic 
scatterers, the ratio of the average power returns from the two scattering classes must be 
maximized. 

In the case of distributed scatterers in the clear, observability is limited only by 
background noise. Since the noise in the receiver channels is statistically independent 
and usually normalized at the same average power, the background noise is generally 
unpolarized. This implies that when the target-to-clutter ratio is maximized, under the 
assumption that the distributed scatterer represents the target, whereas background noise 
denotes clutter, it will be found that the background noise contribution has essentially a 
constant expected power level Giuli. 1985 . Thus, for targets which are assumed to be 
in the clear, or independent of a clutter background, maximizing this ratio is equivalent 
to maximizing the target return power only. Target detection improvement which can be 
attained by such a procedure has been analyzed by resorting to a Gaussian target model 
Giuli. 1982: derived from Huynen’s target decomposition theorem [ Huynen . 1970. 19781. 
Moreover. \an Zyl ct ai. 1987 have determined the optimal co-polarization state for 
maximum power return from an isolated, distributed scatterer. represented by its average 
Mueller matrix. 

For the case of two scattering classes, both of which were either fully or partially 
polarized, Ioannidis and Hammers [1979] employed a Lagrange multiplier method to de- 
termine the transmitting and receiving polarization state that maximized a target’s return 
in the presence of clutter. The target-to-clutter ratio, expressed in terms of average Mueller 
matrices, was maximized in order to determine the optimal transmitting and receiving an- 
tenna Stokes polarization vectors. \an Zyl and Zcbker et al., 1985; 1987] have introduced 
the polarization signature as a means of displaying polarimetric characteristics of various 
scatterers. They Zebker et ai. 1987] numerically determined the polarization state which 
maximized contrast between two classes when the receiving polarization was fixed with 
respect to that of the transmitter, e.g.. co- or cross-polarized returns. 
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It should be noted that all of the techniques previously discussed for polarimetric 
enhancement of a target's return in the presence of clutter have maximized the target-to- 
clutter ratio. However, maximization of the contrast between two classes is not necessarily 
the same as maximizing the targel-to-clutter ratio. For example, the object scattering the 
most power generally is denoted to signify the target class, whereas the other scattering 
class is referred to as clutter, although this may not always be the case. If a target in 
severe clutter is considered, the clutter class may actually scatter more power than the 
target for some transmitting and receiving polarization states. Moreover, the classes can 
exhibit different polarimetric correlation coefficients between the receiver channels while 
having comparable radar cross sections. Thus, the notion of a target and clutter class, 
in some instances, is not well defined. Note also that, in general, maximization of the 
target-to-clutter ratio does not provide the same contrast between classes as maximization 
of the clutter-to-target ratio. Therefore, the problem addressed here is to select the larger 
of these two values and determine its corresponding polarization state. For this reason, 
the term contrast ratio will be adopted as opposed to using either the target-to-clutter or 
clutter-to-target ratio. 

Consequently, the procedure implemented in this paper will determine the transmit- 
ting and receiving polarization state which produces maximum contrast, or separation in 
the average intensity, between the two scattering cksses. To realize this objective, the 
contrast ratio will be maximized, i.e., the maximum contrast ratio is computed in order to 
obtain the optimal linear weighting vector or optimal polarimetric matched filter [ Cadzow . 
19X0: Sovak ct ai, 1987; Swart: 1988 . Processing polarimetric synthetic aperture radar 
(SAR) images with this filter performs a polarization synthesis on the data which yields 
maximum contrast between classes. 

In Section II. the polarimetric matched filter and the contrast ratio are defined. It is 
then shown how to realize the polarimetric matched filter in terms of an equivalent trans- 
mitting and receiving polarization pair. The contrast ratio is defined as a function of a 
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linear weighting vector (polarization filter) and the polarimetric covariance matrices of the 
two classes. The method for achieving maximum contrast between classes, i.e.. determin- 
ing the maximum contrast ratio and is corresponding optimal polarimetric matched filter, 
as well as a closed form solution to this problem for the case in which the polarimetric 
covariance matrices contain four zero elements is considered in Section III. In Section IV. 
the case where the radar transmitting polarization state is fixed will be addressed. Here, 
the receiving polarization state is optimized so that the maximum contrast ratio is at- 
tained under this constraint. Discussion of the results obtained using optimal polarimetric 
matched filtering is the scope of Section Y. 


II. The Polarimetric Matched Filter and Contrast Ratio 


Assume that two classes of statistically distributed scattering types exist. Each class 


.Y 


where j = a.b repre- 


is represented by a covariance matrix of the form Ej = E 
sents class .4 and class B scatterers, respectively. Here E • denotes the expected value 
and superscript f. the complex conjugate transpose operation. For the case of electromag- 
netic waves which are backscattered from a reciprocal media. HV = V H. Therefore, the 
polarimetric feature vector. X. is expressed in a horizontal-vertical polarization basis as 
Kong ct al., 1988 j 


A' = 


HH 

HV 

VV 


( 1 ) 


The objective is to find the best linear weighting vector or polarimetric matched filter 
for processing an observed polarimetric feature vector: that is. the linear combination 


v = n rt a 


( 2 ) 


where 


tr - 


h’aa 

nVr 


(2a) 
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which provides the maximum contrast ratio, r. between the two respective classes (class 
.4 and class B). The maximum contrast ratio is defined as 


r - 


MAx{MAX r (r aJ ).MAX r »r t0 ; 


(3fl | 


where 


r ab - i = 

H -6 »* 


r ba ~ 


ir 

ir f W 


[ 36) 


(3c) 


Note that from a physical point of view, the elements of the vector IV in (2) are linear 
weighting coefficients which adjust the amplitude and phase of the polarimetric radar 
measurements. In equation (3b). the term r a j denotes the contrast ratio of class A with 
respect to class B scatterers. The contrast ratio of class B with respect to class .4 scattering 
elements is expressed in (3c) as r^ a . The symbol MAX{ • } signifies the maximum value 
of the argument, i.e.. either or r^ a in this case, whereas MAX^y(r a £) indicates that a 
linear weighting vector 11' has been obtained which maximizes r a j independently of r^ a . 
and MAX^r(rj a ) implies that a different linear weighting vector has been found which 
maximizes r^ a independently of r a j. Also, the numerator and denominator in (3b) and 
(3c) are obtained from (2) by taking the expected value of the return power from each 
class. 

To demonstrate that the polarimetric matched filter. U\ directly corresponds to 
specific transmitting and receiving polarizations, we express equation (2) in terms of a 
monostatic reciprocal scattering matrix as 


V = Hr Vr 


HH HV 1 


i 

HV IT! 


v t J 


(4) 


in which the values H( and V t . H r and V r represent the horizontal and vertical vector com- 
ponents of the transmitting and receiving polarization state, respectively. Also, without 
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loss of generality, it is assumed that 



H t 2 - \yr = i 

( oa ) 


Hr 2 - Vr 2 = 1 

(56) 

Equating (2) and (4) yields 

-d 

f -c 

*-* 

II 

scT 

(6a) 


H t V r - V t H r = W* hv 

(66) 


VfVr = W’ v 

(6c) 

Thus, given a general linear 
components, Hf. l/. H r . and 

weighting vector, II*. its corresponding polarization 
IV, can be completely specified through (6). This w 

state 
ill be 


shown in Section III. 

Note that equation (4) indicates the transmitting and receiving polarization vectors 
are reciprocal, i.e.. the terms Hf and \) can be interchanged with H r and l' r without 
altering the measurement. Therefore, in the case of reciprocal backscattering. the same 
contrast ratio will be obtained if the transmitting and receiving polarization vectors are 
exchanged. 

III. The Optimal Polarimetric Matched Filter Required to Obtain Maximum Contrast 

Between Two Scattering Classes 

In this section the optimal polarimetric matched filter which produces maximum 
contrast between two scattering classes is determined. It will then be transformed into the 
specific transmitting and receiving polarization state that a radar can utilize to realize the 
same maximum contrast ratio. 

In order to compute the optimal polarimetric matched filter, (3b) and (3c) must 
be maximized. The linear weighting vector which corresponds to the maximum contrast 
ratio, shown in (3a). will be denoted as the optimal polarimetric matched filter. The 
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maximization procedure makes use of the Lagrange multiplier technique. Details of this 
procedure were outlined by Cadzow 1980 . although the steps will be repeated here for 
completeness. For example, in order to maximize r Q in (3b). 

MAX {lT f fair} (To) 

is determined under an arbitrary constraint 

lT t SjiT=l (76) 


This reformulation is possible, without loss of generality, since the linear weighting vector 
can be multiplied bv any arbitrary complex constant without affecting the contrast ratio. 
The solution to this constrained maximization problem is obtained by making use of the 
Lagrange multiplier concept, which reflects the constraint shown in (7b). Its solution will 
be a stationary point of the auxiliary functional 


f(W) = ir f E a H'-A 


1 - \\ A L h U r l 


( 8 ) 


in which A is a scalar valued Lagrange multiplier. Specifically, the stationary points of 
this auxiliary functional are found first by representing the generally complex vector U in 
terms of its real and imaginary components, as - i \V j. Then, taking the gradient of 
the auxiliary functional with respect to and W [, setting the resulting equations equal 
to zero. i.e.. 


*m - o 

mi's 


(9a) 


- o 

air, 


(96) 


yields the necessary condition for a maxima or minima to occur. Carrying out (9a) and 
(9b) leads to the eigenvalue equation 
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Note that the eigenvalue (Lagrange multiplier) A. in (10). is the contrast ratio r a f ) given 
by (3b), whereas 1/A signifies t he contrast ratio r f )a shown in (3c). Since the objective 
is to determine the maximum contrast ratio between classes, the values of the maximum 
and the reciprocal of the minimum eigenvalue must be compared and the larger of the two 
selected. The eigenvector which corresponds to this maxima is the optimal polarimetric 
matched filter that should be employed to process the radar polarimetry. Note that it 
is not required to maximize both (3b) and (3c). By extremizing (3b) then selecting the 
largest of either the maximum eigenvalue or the reciprocal of the minimum eigenvalue, 
both (3b) and (3c) have been simultaneously maximized. 

In the event that the eigenvalues of ( 10) are degenerate, there will exist no preferred 
polarization basis for which the two objects expected power return can be separated. As- 
suming then, that the contrast ratios are not degenerate, the optimal polarimetric matched 
filter. IT. is interpreted to be the equivalent transmitting and receiving polarization state 
which a radar can utilize in order to detect the maximum contrast, or separation in average 
intensity, between classes. 

The contrast optimization approach used for the case of a monostatic radar also can 
be applied to a polarimetric bistatic radar. Taking into account the fact that for bistatic 
scattering HY = \ H when defining X and U in (1) and (2). leads to 4x4 polarimetric 
covariance matrices which characterize the scattering classes. Then applying exactly the 
same method of solution yields the transmitting and receiving polarization state that 
maximizes contrast between scattering classes. 

Once the optimal polarimetric matched filter is obtained, the corresponding trans- 
mitting and receiving polarization state can be calculated. Without loss of generality, the 
case when 11 '/j/ ( is not equal to zero will be shown. From (•">) and (b), it is found that 


H f = 


211 


- 1 


hh 


Ula) 



10 


V'f = H, ■ 


= v ,, r i -’ 2 " 41,1 


116) 


Hr - W 11 66 
< ‘ 


r P = 


H t 


' V hv = \t"to) 2 -M'\h' V rv)' 


where 


1 

v r l • r 2 • r 3 


(llr) 


(11 d) 


( 12 ) 


r l = 



(12a) 


TO - 1 — 


2 "'hh 


U /u - V (U ^ )2 - 4(n '^ n -)‘ 
73=1 ^ 


( 126 ) 


(12c) 


The absolute intensity value. /. is given by 

/ = r lf t T| 2 


( 13 ) 


where r is the amplitude normalization constant given in (12). 

The observed sign change in (11). i.e., - or -. indicates the reciprocity of the trans- 
mitting and receiving polarization state, as previously mentioned. Also, the resulting 
transmitting and receiving polarization state is independent of any multiplicative con- 
stant effecting the matched filter. This is necessary since the general complex eigenvector 
solution to (10) can vary by a multiplicative complex constant; however, the resulting 
polarization state remains unaffected since this constant can be factored out. 
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Finally, a comparison between the methods for contrast enhancement presented in 
this paper versus that originally proposed by loan nidi a and Hammers .19791 is in order. In 
their method, the target-to-clutter ratio was maximized to determine the optimal trans- 
mitting, and receiving antenna Stokes polarization vectors. The loannidis and Hammers' 


method requires the use of three constraints in order to solve the maximization problem. 
One is similar to (7b) in that it constrains the denominator of target-to-clutter ratio to 


be equal to an arbitrary constant. The other two constrain the transmitting and receiving 
vectors to be. antenna Stokes polarization vectors. This results in complex expressions 
which specify the optimal transmitting and receiving polarization state. In addition, they 
do not obtain the matched filter which corresponds to the optimal transmitting and re- 
ceiving polarization state. Using Cadzow's method, only one constraint (7b) is needed to 
solve for the optimal polarimetric matched filter. In this case, the maximization proce- 
dure only requires solving the eigenvalue problem shown in (10). It should be pointed out 
that both methods yield identical results when polarimetric target and clutter classes are 


pre-specified. 

The major difference between these two techniques is that loannidis and Hammers' 
method dealt with the specific problem of maximizing the target-to-clutter ratio by deter- 
mining the optimal transmitting and receiving antenna Stokes polarization vectors, whereas 
the matched filtering approach used in this paper can be applied to a more general class of 
problems. That is. Cadzow's procedure extends to multi-channel, multi-frequency sensor 
data. Polarimetric contrast, enhancement is considered here as a special case. 


Closed Form Solution for the Case of a Covariance Matrix with Four Zero Elements 

Thus far. the most general form of the polarimetric covariance matrix has been 



assumed, which is 
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v _ 


a J 


J j v <j 


pi 


J V ‘J 


J ;vu 


c . 


v *J : J 


— - c a 


L rjy, j 


where j = a. b represents the class .4 ami class B parameters, respectively, and 


12 


( 14a ) 


= *66 
e = <r hv <7 


- c r vr a 


E\HH • IT' 


P = 


rr 


V* 


P exp 



J = 


E\HH ■ HV’ 


<T\j e 


- .jlexp 


‘Oj\ 


C 

s 


\HY • VV* 


v x /c- 


exp 



(146) 
(14c) 
(14 d) 

(14c) 

(14/) 
( 14<?) 


Here, the values <7/,/,. <r/, r . and <7 rr denote the normalized backscatter cross section per 
unit area of the H H. HY . and IT returns Kong e t al., 19N* . 

It has been rigorously shown using the random medium model [ Shin et al., 1986; 
Boryeaud et al.. 19*7’. that when each of the two scattering classes can be modeled as 
a uniform terrain cover, no average correlation exists between HH and HV returns, or 
between VV and HV returns. Therefore, the variables J and £, in (14), are both equal 
to zero and the polarimetric covariance matrices contain four zero elements. This implies 
that the terrain exhibits azimuthal symmetry from a statistical point of view. It should 
be pointed out that this effect has been experimentally observed at. various sites by MIT 
Lincoln Laboratory in their polarimetric measurements at 3o GHz Boryeaud et al.. 1987 
(see Section Y and Table 2). In this case the polarimetric covariance matrix can be 
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expressed as 


_ 't/ 


1 

0 




0 OVV 
0 

0 




( l"») 


A closed form solution to the general eigenvalue problem in ( 10). based on the covari- 
ance matrix in (15), will be now presented. First, the eigenvalues for the matrix Ej, 


are determined. They can be expressed as 


A i — 


n a b 




A o ,3 — 0 

-"bib ( 1 - Pb 2 ) 


la ~ ~b ~ - \'‘b~a Pa\ Pb cos (op. - Op,.) 

la f lb ~ 2 y/lbla Pa\\Pb cos (Op t — Op,.) 
^ - -ba7i> (l ~ iP6i 2 ) ( l ~ Pa 2 ) 


Their corresponding eigenvectors are given as 

Ai-ir 1 = 

A 2 .3 — H' 2.3 = 


1 

0 

*2.3 


where 


T _ l i 

t 2.3 “ 


la - lb - ' - v'lbla Pa Pb\ sin (op.. ~ Op, ) 

[*6 ~ *a “ - v 'TfcV a Pa Pb cos (Op., - Op,.) 
\ - 4* a* 6 l 1 ~ Pb ') ( 1 ~ Pa 2 ) 


(ibv * aPa - la v 


(16a) 


(166) 


(17a) 


176) 


(17c) 


As in the case of the generalized solution, the eigenvector corresponding to the max- 
imum of either the largest or the reciprocal of the smallest eigenvalue in (16) will produce 
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maximum contrast between the two classes. Therefore, this eigenvector will be the opti- 
mal polarimetric matched filter which should be utilized to process the polarimetric feature 
vector. 


IV. The Optimal Receiving Polarization State for a Fixed Transmitting Polarization 

If presented with a situation where the radar transmitting polarization state is fixed, 
an optimal receiving polarization state can be determined which maximizes the contrast 
between the two classes. This problem arises in the case of radar systems which are not 
fully polarimetric. i.e.. they transmit using only a single polarization, say horizontal, but 
receive the principle and cross-polarization components of the scattered response, say the 
horizontal and vertical returns. Applying the following technique will indicate how to 
coherently combine the horizontal and vertical returns such that the contrast between 
classes is maximized. 

Assume that the values Hf and \) are known and that the requirement shown in 
(5a) is satisfied. Thus, from (6), IV may be written as 


IT = 


Substituting (18) for U' in (3b) gives 


(H t H r )* 


\ H t 

o ' 


Hr' 
. V r . 

(H t V r ) m -(V t H r )* 

— 

If" 

Ri 


w>r 


L 0 



= < R 


( 18 ) 


R* Z a R 

r ab = 

V Zb R 


( 19 «) 


where 


Z a = C 

= C f Sft c 


(196) 

(19c) 


Since Z n and Zf, are hermitian symmetric, positive semidefinite matrices. (19a) can be 
extremized as in the previous section, to obtain the generalized eigenvalue problem 


Z a R = \ 'Z b R 
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Once R is obtained, it should be normalized so that (5b) holds. Application of the trans- 
form given in (18) will yield H . 


V. Results and Discussion 


In order to present the optimal polarimetric matched filtering results in a compact 
format, orientation U') and ellipticity (\) angles (Fig. 1) are utilized to express the 
transmitting and receiving polarization states. Here, the definitions from Kong 1986i are 
adopted. Therefore, horizontal ( H ) and vertical (\) polarization states will have zero 
degree ellipticity angles, with orientation angles of O'" and 90°. Right (R) and left {L) 
polarization states are orientation independent with ellipticity angles of 45° and —45°, 
respectively. In addition. 0" < v < 1*0° and —45° < \ < 45 r . 

A general polarization vector. 


{Ph} 


\ p h\ exp [io h ] 

[Pv\ 


_ exp [io v | _ 


written in terms of horizontal ( h ) and vertical (v) vector components, can be transformed 
into a normalized Stokes vector as follows 


' S 0 ' 


' Ph ' 2 - A ! 2 ' 


1 

- S 'l 


Ph - - Pv i 2 

2 P^\\P V \ coso 

= ( Ph 2 * Pv 2 ) 

cos 2\ cos ' 111 ' 
cos 2\ sin 2c- 

•** 3 J 


. 2 P/, Pv sin<t>. 


sin 2\ 


where o - o v - o I’sing this equation, the angles c and \ are obtained. 

Two data bases were utilized to study the contrast problem, utilizing the techniques 
outlined in the previous sections. Table 1 gives polarimetric covariance statistics extracted 
from San Francisco Bay area. L-band (1.225 GHz) SAR data, collected by the Jet Propul- 
sion Laboratory's airborne polarimeter [Ztbk(re\ al., 19X7 . These covariance statistics 
were obtained from the training areas shown in Figure 2 and were utilized to generate the 
results shown in Tables 4. 5. and 7. Similarly, the experimental polarimetric covariance 
data, shown in Table 2. was supplied by the MIT Lincoln Laboratory Borgcaud et al., 
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19S 7: Kong c t ai. 19S<S _. This data, collected at 35 GHz. was used to generate the results 
presented in Table 6. The MIT Lincoln Laboratory radar imaged a vegetation field consist- 
ing of grass or trees at a range of approximately 2 kilometers. Studies using this database 
indicate that essentially no correlation exists between the HH and H \ . and between the 
HY and 11" polarimetric returns, i.e.. the terrain clutter exhibits azimuthal symmetry. 

Next, a comparison is made which evaluates the performance obtained by processing 
polarimetric radar data using the optimal polarimetric matched filter versus other com- 
monly used polarization states. Linear weighting vectors which correspond to commonly 
used transmitting and receiving polarization states are given in Table 3. The weighting 
vectors presented in this table were generated using equations (6a) through (6c). (21) and 
(22). These linear weighting vectors are expressed in a horizontal- vertical polarization 
basis. 

Table 4 presents theoretical contrast ratios r a ^ and rj )a obtained when utilizing the 
above mentioned transmitting and receiving polarization states (Table 3) as well as the 
optimal solution. Here, class .4 and B scatterers have been defined to denote the park 
and urban (city) regions, which were represented by their corresponding covariance ma- 
trices. Ha and H/,. respectively. As previously discussed, for reciprocal backseat tering the 
transmitting and receiving polarization state may be interchanged while maintaining the 
same contrast ratio. This is clearly indicated by the HY and Y H results and the LR 
and RL results. The values r ^ and r^ 0 . which denote the maximum and reciprocal of 
the minimum eigenvalues found after solving (10). are expressed in terms of their corre- 
sponding orientation and ellipticity angles, e and y. From Table 4 it is seen that the 
maximum contrast ratio between the two selected classes is 9.3X dB. Note that had only 
r n i been maximized, a contrast ratio of 2.37 dB would have been realized. In some cases, 
though, this may be what is required. If the problem was only to make the park processed 
pixel intensity as large ss possible with respect to that of the city, the transmitting and 
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receiving polarization state corresponding to the average power ratio of 2.37 dB would be 
the appropriate matched filter which should be used to process data. 

Table 5 shows the actual polariniet ric contrast enhancement achieved when process- 
ing the radar measurements using the optimal polarimetric matched filter and various 
other polarization filters. Since the contrast ratios for r^ a . given in Table 4, are larger 
than those for r a £, the linear weighting vectors which correspond to the contrast ratios for 
r^ a have been used to generate the results shown in Table 5. Thus, Table 5 contains the 
actual average processed pixel intensity realized for each of the two classes (urban and park 
areas) for both suboptimal polarization filters, i.e.. transmitting and receiving polarization 
states which do not provide maximum contrast between classes, and the optimal polari- 
metric matched filter. In comparing the data in this table, it is seen that the quantitative 
measure of attainable contrast is the contrast ratio, which is the linear ratio of (or the 
logarithmic distance between) the average pixel intensity for the two respective classes. 
In the case of processing data with the optimal polarimetric matched filter, this distance 
is maximum. Thus after optimal processing, it is possible to more readily separate the 
two classes than prior to it. since the distance between the average value of pixel intensity 
has increased between the two classes. Also note that the contrast ratios, shown in Table 
5, and those given in Table 4 for r^ a . are similar indicating a good match between the 
theoretical predictions and the processed results. 

A further demonstration of contrast enhancement can be seen visually by comparing 
Figures 3A through 3D which show San Francisco Bay area images. This imagery has been 
synthesized utilizing some commonly employed linear polarization states, in addition to 
the transmitting and receiving polarization state required to produce maximum contrast 
between the park and urban regions. In all four of these images, the average processed pixel 
intensity of the park region was set to the baseline value, i.e, the minimum quantization 
level of the imaging system display, which was -20 dB. The maximum quantized intensity 
was -10 dB. By utilizing these quantization limits, some clipping of the higher and lower 
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intensity levels has occurred. This procedure was implemented to compare more easily the 
contrast between images. Figure 3A illustrates the result of processing the San Francisco 
Bay polarimetry using the optimal polarimetric matched filter. The contrast ratio obtained 
between the city and park area, as previously indicated in Table 5. was 9.4 dB. Contrast 
ratios achieved using HH (Fig. 3B). l*V ( Fig. 3C ). and HY ( Fig. 3D) polarization filters 
were only 7.3. 5.4. and 2.6 dB. respectively. As indicated in Table 5. utilizing the optimal 
transmitting and receiving polarization state, i.e.. the optimal polarimetric matched filter, 
to process data yields the maximum contrast ratio. Figure 3 also shows that the optimal 
and H H synthesized images appear somewhat similar; this is due to the fact that there is 
only a 2.1 dB difference between their contrast ratios. However, the optimal polarimetric 
matched filter always yields a larger contrast ratio between classes than when any other 
transmitting and receiving polarization states are utilized. 

Contrast ratio results, obtained using the MIT Lincoln Laboratory data are presented 
in Table 6. As was the case in Table 4. theoretical contrast ratios are given for frequently 
employed polarization states as well as for the optimal solution. In this table, the tree and 
grass regions have been arbitrarily selected to denote class a and 6 scatterers. respectively. 
The optimal solution again is represented by its corresponding orientation and ellipticity 
angles v and \. in which case, the values presented for r a ^ and r^ a . shown in Table 6, 
signify the maximum and reciprocal of the minimum eigenvalues found when employing 
the equations shown in (16a) and (16b). 

In Table 7, the optimal receiving polarization states required to produce the maxi- 
mum contrast ratio between classes for various fixed transmitting polarization states are 
presented. These results show that by employing the optimal receiving polarization state, 
all contrast ratios have increased relative to those shown in Table 4 or 5. Thus, for a given 
transmitting polarization state, synthesizing imagery using the optimal receiving polar- 
ization state always yields a larger contrast ratio than when any other polarization state 
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is used. However, the maximum contrast ratio is achieved only when both the optimal 
transmitting and receiving polarization state is employed. 
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Figure 2 


Figure 3 


List of Figures 


Generalized elliptic polarization state. 

Training areas used to generate the covariance matrix elements for the park 
and urban (city) regions, shown in Table 1. 

San Francisco Bay area images synthesized using t lie optimal polarimetric 
matched filter (A I. H H ( B h W (C ). and HY ( D ) polarization filters. The 
corresponding contrast ratios between the city and park region were 9.4. 7.3. 
o.4, and 2.6 dB. respectively. 
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rr(dB) 

1 6 

j 1 P 

-J 

Op 

3 

°3 


Oc 

s 

Park 

-40.5 

0.406 

j 1.42 

0.219 

-0.463 

0.168 

0.590 

0.090 

-1.22 

Urban 

-33.2 

0.137 

'< 0.907 

0.404 

2.98 

0.792 

0.001 

0.414 

-2.98 


Table 1 Covariance 

matrix elements for 

park and 

urban (city) 

regions. 

Phase angles 



are given in radians. 

rr\dB) 

t' 

-J 

P 

Op 

Trees - 13.0 

0.06 

1.1 

0.74 

0.0 

Grass - lo.Q 

0.15 

1.2 

0.56 

0.0 


Table 2 Covariance matrix elements for a uniform terrain cover consisting of grass and 


tree regions. 
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Transmitting Receiving Linear Weighting 


Polarization Polarization Vector : IT 


H H 

'1.0' 

0.0 

0.0 

' 

H V 

'o.o' 

1.0 

0.0 

I 

r h 

'o.o' 

1.0 

0.0 


r v 

'o.o' 

0.0 

1.0 


L L 

0.5 

i 

-0.5 


L R 

'O.o' 

0.0 

_0.o_ 


R L 

'O.o' 

0.0 

O.o 

1 

R R 

O.o 

-i 

-O.o 



Table 3 Commonly utilized transmitting and receiving polarization states versus their 
corresponding linear weighting vectors (expressed in a horizontal-vertical po- 
larization basis ). 
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Transmitting 

Receiving; 

Contrast Ratio (dB) 

Contrast Ratio (dB) 

Polarization 

; Polarization 

r ab 

r ba 

H 

H 

-7.30 

7.30 

H 

r 

-2.58 

2.58 

V 

H 

-2.58 

2.58 

V 

V 

-5.35 

5.35 

L 

L 

-6.94 

6.94 

L 

R 

-3.29 

3.29 

R 

L 

-3.29 

3.29 

R 

R 

-6.73 

6.73 

v = 1.82 

v = 107.0 

2.37 

-2.37 

X = 3.72 

\ = -1.64 


1 

r = 48.7 

c = 150.3 

-9.38 

i 

! 

9.38 ; 

V = -6.44 

X — 3.51 




Table 4 Theoretical contrast ratios between classes when U. a = Park and Ej, = Urban 
San Francisco Bay regions, respectively. Orientation and ellipticity angles are 
given in degrees. MAXj^ (r a j) = 2.37 (IB. MAX|pr(r^ a ) = 9.38 dB, therefore 
r = 9.38 dB. 
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Transmitting 

Polarization 

Receiving 

Polarization 

Urban 

Class 

(dB) 

Park 

Class 

(dB) 

Contrast 

Ratio 

(dB) 

H 

H 

-33.2 

-40.5 

7.3 

H 

V 

-41.8 

-44.4 

2.6 

V 

V 

-33.6 

-39.0 

5.4 

L 

L 

-34.2 

-41.0 

6.8 

L 

R 

-38.6 

-41.9 

3.3 

R 

R 

-34.0 

-40.8 

6.8 j 

v — 18.7 

c - 160.3 



! ] 



-34.3 

-43.7 

9.4 ; 

\ = -6.44 

X = 3.61 



i ! 


Table 5 Actual average pixel intensities and the contrast ratios between the park and 
urban San Francisco Bay regions when data was processed using commonly 
employed polarization filters and the optimal polarimetric matched filter. 
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Transmitting 

Receiving 

Contrast Ratio 

Contrast Ratio 

Polarization 

Polarization 

r ab 

r ba 

H 

H 

2.00 

-2.00 

H 

V 

-1.9s 

1.98 

r 

H 

-1.98 

1.98 

V 

V 

1.62 

-1.62 

L 

L 

-1.00 

1.00 

: l 

R 

2.28 

1 -2.28 

R 

L 

2.28 

i -2.28 

R 

R 

-1.00 

1.00 

c — 0 

c = 90 

-1.98 

1.98 

\ =0 

\ =o 



— 0 

c ^ 0 

2 . n 

i i 

! -2.31 

y = -38.3 

\ - 38.3 


: 


Table 6 Theoretical contrast ratios between classes when trees = a and grass = b for 
class .4 and B scatterers. respectively. Orientation and ellipticity angles are 


given in degrees. MAX^?(r a ^) = 2.30 dB. MAX jjr ( rj, a ) = 1.99 dB. therefore 
r = 2.30 dB. 
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Table 7 
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Transmitting 

Receiving 

Contrast Ratio 

Polarization 

Polarization 

(dB) 


c = 31.8 


H 

\ -8.04 

7.83 


t- - 134.2 


V 

\ = 4.34 

6.06 


r = 27.0 


R 

\ = 26.1 

6.97 


f = 169.1 


L 

X = -21.4 

7.36 

v = 48.7 

v = 150.3 

i 

9.38 

1 

X. = “6.44 | 

X = 3.51 



Optimal receiving polarization state for a fixed transmitting polarization state 
when H a = Park and H/, = Urban San Francisco Bay regions. Orientation and 


elli pticity angles are given in degrees. 



m 

v 



Figure 1 Generalized elliptic polarization state. 
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ure 2 


Training areas used to generate the covariance matrix elements for the 
park and urban (city) regions, shown in table 1. 
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Figure 3 San Francisco Bay area images synthesized using the optimal polarimet- 
ric matched filter (A), HH (B), VV (C), and HV (D) polarization 
filters. The corresponding contrast ratios between the city and park 
region were 9.4 , 7.3 , 5.4 , and 2.6 dB, respectively. 



